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Abstract 



(N 
> 

o 
o, 

\^ , We use Coulomb gas methods to derive an explicit form for the scaling limit of 

^^ ' the partition function of the critical 0(n) model on an annulus, with free boundary 

conditions, as a function of its modulus. This correctly takes into account the 
magnetic charge asymmetry and the decoupling of the null states. It agrees with 

1 -Q I an earlier conjecture based on Bethe ansatz and quantum group symmetry, and 

C^ ' with all known results for special values of n. It gives new formulae for percolation 

(the probability that a cluster connects the two opposite boundaries) and for self- 
avoiding loops (the partition function for a single loop wrapping non-trivially around 

1^ - the annulus.) The limit n — > also gives explicit examples of partition functions in 

\^ . logarithmic conformal field theory. 



*Address for correspondence. 



1 Introduction. 

The Coulomb gas approach to two-dimensional critical models which can be written as 
loop gases, such as the 0(n) model and Q-state Potts models, has been extraordinarily 
successful in deriving their universal bulk properties. First developed by den Nijs Ij and 
Nienhuis in order to explain conjectured exact values for the principal bulk critical 
exponents, it was adapted by di Francesco, Saleur and Zuber |^ to compute the partition 
function on the torus, which encodes all the bulk scaling dimensions jl|. 

However, this particular approach has not, so far, been successfully adapted to explain the 
conjectured exact values E] for the boundary scaling dimensions in these models, even 
though many of them have now been derived using other methods "H] applied to related 
lattice models, and, more recently, rigorously using Schramm-Loewner Evolution (SLE) 
[7j. From the point of view of conformal field theory (CFT), both the boundary and the 
bulk exponents are encoded in the partition function on the annulus jH], and therefore 
a direct computation of this object is of great interest. From the point of view of CFT, 
general values of n and Q give irrational non-unitary examples. 

As we shall explain, the naive application of Coulomb gas methods to domains with 
boundaries fails to account either for the reduction of the central charge c from its free 
field value of unity, or for the elimination of the null states, which, as is known from CFT, 
is necessary to maintain modular invariance. In this paper we propose a resolution of 
these issues which provides an explicit formula for the partition function on the annulus. 

Consider an annulus {0 < x < i,0 < y < L), identifying x = 0,i. Introduce the conjugate 
moduli 

Note that if the annulus is conformally mapped to the region Ri < r < R2 between two 
circles, q = ln(i?2/-Ri). We impose free boundary conditions on the 0(n) spins on y = 0, L. 
As usual in the Coulomb gas, we introduce the parametrisation n = \fQ = 2 cos x, 
g = 1 — x/tt, where 1 < g < 2 corresponds to the dilute critical point of the 0(n) model 
(or the tricritical point of the Q-state Potts model), and | < (7 < 1 to the critical dense 
phase of the 0(n) model (or the ordinary critical point of the Potts model.) Then our 
main result for the annulus partition function is 

-^TT/1 rx-i Y^ sin(p-M)x aZ 
^ -*■ ^-^ sm X 



r=l 



In terms of the conjugate modulus q this becomes 

/ , Ni/9 <= TT/ 9^N 1 v^ sin((Y + 2m7r)/gf) (x+2^mf {i~gf 
Z={2/gy/'q--^l[{l-fT'Y.^^^- '-^q^^^ '^ . (2) 
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This has the form expected from boundary conformal field theory jH] : (pj is 

Z = q~'^^^'^'^/^dAXAiq), where A runs over the allowed set of boundary scaling dimen- 



sions, Xa{<i) is a highest weight Virasoro character, and dA is a degeneracy factor, which 
is a polynomial in n (although only integer in those cases when the theory is unitary). In 
this form the explicit expression (^ is not new, and indeed is originally due to Saleur and 
Bauer Pj, who deduced the allowed scaling dimensions from Bethe ansatz and quantum 
group arguments. In this context the degeneracy factor d/s. is the quantum dimension. It 
has been used extensively in papers by Saleur and Pasquier jTU], Saleur [TT], and more 
recently appeared in a paper by Read and Saleur ^j . However, a direct derivation from 
the lattice 0(n) model solely using Coulomb gas arguments has not appeared to our 
knowledge, and this is the point of the present paper. 

Note that Q also has the form expected from boundary CFT [9^ , namely 

~-c/i2 ^^ I&aPxaI'?^) where now the sum is over allowed bulk scaling dimensions 2A, 

and 6a is a matrix element with a boundary state. In particular, for A = we have 

^.^_^ygy,.MEM^ (3) 

SmTT^f 

which gives the boundary entropy T^ Infco- 

We have explicitly written the dependence on x i^i Eqs. fll|2|) because, if we wish to 
consider a modified partition function in which the loops which wrap non-trivially around 
the annulus are counted with a different weight n' = 2cosx', it is simply necessary to 
replace x ^ x'- This allows us to compute interesting quantities for percolation and 
self-avoiding loops. 

For example, in critical percolation {Q = 1) the probability that a cluster connects the 
two boundaries of the annulus is 

oo 

^ = n(i-^r^E(^'^-'^-^'^"""^)- (4) 

r=l fceZ 

The partition function for a single self-avoiding loop which wraps non-trivially around 
the annulus (the number of such loops weighted by fi^^^'^^^^, where /i is the non- universal 
connective constant) is 

oo 

Zi = ll(^-qT'j:H-ir'^'^''-''^- (5) 

r=l fcez 

In the limit g — >■ 0, we find Zi ~ (l/67r)| lng|. The form of this agrees with a rigorous 
result of Werner ^3] . 

The layout of this paper is as follows. In Sec. |21 we give a brief survey of Coulomb gas 
methods as applied in the plane and the cylinder, and then discuss the particular problems 
associated with domains with boundaries, in particular the annulus. This will lead to the 
proposal (HJ for the partition function. As with most Coulomb gas methods, this is not 
wholly deductive, but rehes on some heuristic reasoning. However, in SecISl we show that 



(^ agrees with previously known results for various special cases (for example the Ising 
model and the 3-state Potts model.) In Sec. HI we then derive a variety of new results, 
some of which have already been mentioned above. 



2 Coulomb gas on the cylinder and the annulus. 

2.1 Basics. 

We first summarise the Coulomb gas arguments as applied to the plane and cylinder, as 
formulated by de Nijs [IJ and Nienhuis 0, and elaborated by Kondev [T3] . 

The 0(77,) model is most easily realised on the honeycomb lattice. At each site r is an 
n-component spin s(r) (initially ra is a positive integer.) The Boltzmann weight for a 
given configuration is 

~(l + ts(r)-s(r')), (6) 



n 

r,r' 

where the product is over all edges (r, r') of the lattice. The partition function is the 
trace over these weights, a linear operation defined by Tr 1 = 1, TrSa(r)sfe(r) = 5ah and 
Tr Sa(r) = Tr Sa(r)sfe(r)sc(r) = 0. Expanding ® in powers of t gives a sum over all 
subsets Q of the edges, with an associated factor t'^L Implementing the trace operation 
eliminates all subgraphs which are not unions of non-intersecting closed loops (for the 
time being we ignore boundaries), and each of these gets counted with a weight n. 

At this point we can allow n to be any positive real number. This gives a measure on 
the allowed subgraphs Q, called the loop gas. If t is small, the mean loop length is finite, 
even in the thermodynamic limit, but there is a critical value tc at which it first diverges. 
This is called the dilute critical point. For t > tc a. single loop contains a finite fraction 
of the sites: this is the dense phase. 

The critical Q-state Potts model on the square lattice can also be written, via the Fortuin- 
Kasteleyn [16j correspondence, in terms of a loop gas, in which each closed loop carries a 
factor \/Q. 

Both these loop gas models can be mapped to a model of heights h{R) on the sites R of 
the dual lattice, by first orienting each loop, so that a configuration of A^ non-oriented 
loops corresponds to 2^ configurations of oriented loops, and then, for each edge of the 
lattice, assigning height differences Ah = 0, ±7r between the neighbouring sites of the 
dual lattice according to whether the edge is contained in the oriented subgraph Q, and 
its orientation. The weight n (or a/Q) for each non-oriented loop is distributed into a 
factor e^*^ for each clockwise (anticlockwise) oriented loop, where n = 2cosx- Although 
these weights are complex (a feature which lies at the heart of the difficulties associated 
with a rigorous treatment of the Coulomb gas approach), they have the advantage of 
being local, in the sense that they may be distributed so that each loop acquires a factor 



g*ex/2T whenever it turns through an angle ^ at a vertex. 

However, it should be noted that, at least for the fully packed model on the square lattice, 
there is a mapping to the 6-vertex model with positive Boltzmann weights: at each vertex 
the two loops are either oriented parallel to each other, with weight e*^'^ ■ e~*^'^ + e~*^' ^ • 
gix/4 = 2, or anti-parallel in which case the weight is 

eix/4 . gix/4 ^ g-ix/4 . g-ix/4 _ 2 cos(x/2) = {n + 2f'^ . 

The Coulomb gas method assumes that, in the continuum limit, the discrete heights 
become continuous and the Boltzmann weights converge to e~^ where S is the action of 
a free field theory 

S = (g/An) [{dhf dxdy . 



The original discrete model may be recovered from this by adding a term A YIr cos 2h{R) 
in the limit A — i> — oo. 

However, on a cylinder of length i and circumference L, with i ^ L, this does not 
properly account for loops which wind around it: these can be taken into account by 
placing 'electric' charges e^*'^^/'^)^ at either end. This modifies the partition function to 
Z ~ e'^'^'^^, identifying the total central charge 

c = 1 — . 
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The scaling dimensions of electric charges e*'''* are also modified if we put them at the 
ends of the cylinder as well: 

Note that Xg ^ x_g: this an example of the electric charge asymmetry introduced by this 
construction. 

g is fixed in terms of x by requiring ^2] that cos2h be marginal in the sense of the 
renormalisation group, i.e. X2 = 2. This fixes 

g = l± (x/tt) , 

with the sign depending on whether we choose X2 or a;_2. In fact this ambiguity is to be 
expected: for each value of n, x is only defined up to a sign (actually we can add multiples 
of 277 as well, but these give less relevant operators) and these correspond to the dilute 
{g > 1) and dense {g < 1) cases of the critical 0(n) model. In the following we take the 
lower sign by convention. 

Note that these ideas are easy to extend to the partition function on the torus, correctly 
taking into account loops which wrap around some combination of the two cycles 0. 



Now consider the case of the annulus. Throughout this paper we assume free boundary 
conditions on the 0(n) spins, which means that there are only closed loops in the loop 
gas representation. (Partial results using Coulomb gas methods were found for the case 
of fixed boundary conditions in Ref. fT?! for the special case ra = 1 in the dense phase.) 

First consider the case when i ^ L, where we expect 

Z ~ e ' ^ q ' , 

with c given as above. In this limit there is no contribution of loops wrapping around the 
annulus, so we expect that that h{y = L) = h{y = 0). Naively then, we get a free field 
theory with (equal) Dirichlet boundary conditions, which gives c = 1. 

Where does the correction to c come from? One (incorrect) possibility is as follows: 
looking back at the lattice construction, we see that there are extra factors of e^*-*^/^ 
whenever a loop is next to the boundary, which are not properly taken into account in 
the the bulk Boltzmann weights of the height model. The sign is determined by whether 
the height at a site next to the boundary is ivr. In the continuum limit these would lead 
to boundary terms in the action proportional to ix J d±h dl where d±_ is along the inward 
pointing normal to the boundary and dl is a line element. For the annulus these would 
give something proportional to 



ix I {dyh{x, y = 0) - dyh{x, y = L))dx (7) 

However, an explicit calculation (see Appendix) shows that such a combination does not 
contribute to c. In fact, if we add to the action a general boundary term 






/ {aidyh{x, y = 0) + a2dyh{x, y = L))d: 

we find that the effective central charge is 

c=l-{2A/g){ai + a2f. 

Thus not only is there no contribution if ai = —0:2, as in ^, we must also have cti + 02 
real, rather than pure imaginary. 

An equivalent, and easier, way of getting the same modification to c is to assume that the 
correct boundary conditions, even when i ^ L, are 

h{y = L)-h{y = 0) = 7imoj^O. 

In that case we can write h = nrnQy/L + h, where h vanishes on both y = and y = L. 
The functional integral over h gives c = 1 as before, and the modification to the partition 
function is ~ exp{—{g/4:7!-){'7imQYi/L). So if we take 

mo = ±x/t^9 , (8) 






Figure 1: The screening of magnetic flux in a long rectangle. If the background flux m is too 
large, vortex pairs of strength ±2 can be shed from either end of the rectangle and will annihilate 
in order to reduce the free energy. If m is too small, the opposite effect occurs. 



we get the correct c. Note that, for a long rectangular strip rather than an annulus, this 
is like adding magnetic charges, or vortices, at the ends. 

Thus it would appear that there should be a spontaneous average magnetic flux around 
the annulus, ii g 7^ 1. This may be understood heuristically in terms of the preferred 
parallel, rather than antiparallel, alignment of neighbouring loops: this effect should be 
enhanced near a boundary, since loops are geometrically constrained to lie approximately 
parallel or anti-parallel to the boundary. 

The actual value of mo may be fixed by the following argument. Consider first the 
geometry of a long rectangle with i ^ L, where the loops are allowed to end on the 
boundaries at a; = 0,i, but, as before, not on y = 0,L. In that case the total charge m 
flowing along the rectangle is not fixed, and we can ask the question what is its mean 
value tuq in the state of lowest free energy. A flux m corresponds to vortices of strengths 
±m at the ends of the rectangle. However this can increase or decrease in units of 2 by 
shedding vortices from either end of the rectangle (see Fig. Q). The additional free energy 
for creating such a pair of vortices is (g/Arr) ((m + 2)^ — m^) (n/LY{£L) scaling dimension 



of a vortex of strength ±2, from which we read off the scahng dimension 

A2 = ((7/4)((m + 2)2-m2) 

If m is too large, A2 > 1, which means that the corresponding renormahsation group 
?/2 = 1 — A2 < 0. (Note that since we are in the boundary, rather than the bulk, situation, 
the eigenvalue is 1 — A rather than 2 — A.) This implies that such ±2 vortex pairs are 
closely bound. Thus any such pairs shed from the boundaries at x = 0,i will annihilate 
to reduce their free energy. On the other hand, if m is too small, 1/2 > 0, which means 
that any such vortex pairs will unbind. This will act to increase the effective value of m. 

This screening effect implies that the mean value ttiq of the magnetic flux which minimises 
the free energy corresponds to 2/2 = 0, that is 

{g/A){{mo + 2Y-ml) = l, 

so that 

mo = {l-g)/g = x/T^9, 

which is the same as found above in (jH)). There is a similar minimum free energy solution 
with m = —mo. Note that the above argument is analogous to the earlier argument which 
fixed g, where we demanded that electric ±2 charges should be marginal. 

If we now go the annulus, we expect a total average magnetic flux ±mo to spontaneously 
form, even when i ^ L. Now suppose that i/L is not so large, so we can have extra loops 
wrapping around the a;-cycle. We can orient these as before. If the total number of up 
arrows minus down arrows (the additional magnetic flux flowing along the annulus) is p, 
then we get h{y = L) — h{y = 0) = 7r(p ± mo). As for the cylinder, in order to count 
them correctly we need to put in a factor exp{i{p±mo)x). Thus we get the following flrst 
guess for the partition function on the annulus: 

Z = Zo^e'(P+'"»)>^e-(3/^)("P+'"°)'("'/^) + (mo ^ -mo) , 

where Zq = Q~^^'^^Y['^ii^ ~ Q^)'^ is the partition function from h. Note that we should 
sum over both possible signs for mo. If we let p —>■ — p in the flrst term this simplifles to 



00 



Z = g-^/24 JJ(l _ g-)-! Y^ cos{{p - mo)x) g^^/^)^'-^!-^)^/^ _ (g^ 
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Note that if we want to count loops wrapping around the annulus with a different weight 
n' = 2cosx', we just change x ^ x' ^^ the above (keeping g the same.) 

Eq. Q has some good features and some bad ones. In general we expect that Z can be 
written as a sum of terms q^ where h runs over all the allowed scaling dimensions of the 
allowed boundary operators. (For a unitary theory the coefficients should be non-negative 
integers, but this doesn't have to hold for general n.) We see in Q for p = A^ > 1 the 



scaling dimensions of the boundary A^-leg operators, as first conjectured by Saleur and 
Duplantier 0. 

However, for the dilute case with g > 1, p = —1 actually gives the next-to- leading 
term as q ^ 0. This doesn't make sense: we expect this to come from p = N = 1. More 
seriously, (jHl) fails to account for the fact that the scaling dimensions of the boundary A^-leg 
operators correspond to those of the degenerate cases /ii,Ar+i of the Kac table: in general 
these operators correspond to highest weight states whose Virasoro representations are 
reducible: they have a null descendent state (which corresponds to a term in the expansion 
of Ylri^ ~ Q^)~^) of dimension /ii^at+i + A^ + 1. For a unitary theory, or more generally a 
minimal model, we know that such states (and all their descendents) should be subtracted 
out of the partition function. For a non-unitary theory this is not necessary: however we 
shall show later that if they are retained, the behaviour as g — > is incorrect. This leads 
to the conclusion that each term in Q should be modified according to 

Note that this has the feature of automatically eliminating the 'rogue' state aX p = —1. 

We now give a physical argument for this subtraction. Once again it is useful to think 
about the rectangle geometry where magnetic flux can be created or destroyed at the 
boundaries at x = 0,i. It is also useful to think in terms of the energy eigenstates of 
the hamiltonian {'k/L){Lq — c/24) which generates translations in x. In general, each 
configuration with total magnetic flux p will be accompanied by excitations of the h field, 
which correspond to the Virasoro descendents. If the energy of these is correct they can 
resonate with the original highest weight state plus a number of pairs of marginally bound 
±2 magnetic charges, each of which has energy vr/L. If the excitation energy is {p+l)Ti / L, 
exactly p + 1 vortex pairs can be shed from the boundaries at x = 0, £ (see Fig. |2I). These 
states, however, are identical to those with total magnetic charge p — 2(p + 1) = — 2 — p, 
and should therefore not be doubly counted. 

The effect of this subtraction is therefore to modify the sum in (^ to 



y^ cos((]9 - mo)x) (g'^i-J'+i - g'^L-p-i) . 



pez 
Relabelling p —>■ —2 — p in the second term has the effect of modifying 

cos{{p - mo)x) -^ cos((p - mo)x) - cos((p + 2 + mo))x) oc sin((p + l)x) , 

which finally leads to the conjecture (P) for the annulus partition function, after normal- 
ising so that the coefficient of the p = term (the contribution of the identity operator) 
is unity: 

-^TT/1 rN-i Y^ sin(p+ l)x sZ 
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Figure 2: Mechanism for the appearance of nuh states. In this example, an excited state in the 
p = 1 sector has just sufficient energy 27r/L to allow two marginally bound pairs of ±2 vortices 
to form. These can then move to the ends of the rectangle, and one of them can then annihilate 
with the original flux line. The state is therefore equivalent to the ground state in the p = —3 
sector, which has already been counted in Q and which therefore must be subtracted off. 

Note that all the coefficients in (^ are polynomials in n = 2 cosx as we expect. For p = 1 
we get exactly n (the degeneracy factor for a single loop wrapping around the annulus), 
for p = 2 we get n^ — 1, and so on. 



2.2 Modular properties 

Now we express (^ (for general x') in terms of the conjugate modulus q = e~'^'^^'^. Setting 
q = e~^ we have 

Z = Zog^(sinx')-'lme^^' ^^^'^'^"'^"^"^^^ • 
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Using the Poisson sum formula ^ — > ^^ / dpe'^'"^"^^, the integral is 



'p 

1 N./ I o_^ .■(l-9)'^"l2 



The last term in the exponential cancels the g(^~'=)/24_ Under a modular transformation 

Zo = {6/2iTf'^~q-^^ \[{l - q' 

so we end up with 



~2r\-l 



z = (2/^)^/^r^ 11(1 - er' E ^^ -" /^^^^ g 



2 , 2 



sin((2mx + x')/9) _i2L+2iiq)l^A. 



sinx' 



2vr^ 



which finally simplifies to (j21). Note that if we had not subtracted off the null states, as 
in Q, we would find Z ~ g^^/^^ rather than g~^/^^. This is an example of how the null 
state subtraction is necessary to maintain the correct modular properties |3]. 

The leading term as g — > 0, with m = 0, agrees with the exponent found in Ref. ^H] for 
the case of loops wrapping around a long cylinder counted with weight n' = 2cosx'- We 
now have also the prefactor: 

z~(2/,)V2!^!M^,~^. (10) 

sinx' 

If we set x' = X; the other exponents in Q are those of even electric charge operators 
X2m- For general g these are not in the Kac table, consistent with the fact that there is 
no explicit substraction of null states in (0), so the characters are simply given by the 
infinite product. 

We now check (^ against some known cases. 



3 Comparison with known results. 

3.1 n = in the dilute regime. 

In this case g = 3/2, x = ~^/2, for which we expect Z = 1, since loops wrapping round 
the annulus should all get a weight n = 0. ((H) gives 



^ = 11(1 -^')"'E^i^((P + 1)^/2)^ 



3£_ , p 
8 """4 
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The prefactor is +1 if p = (mod 4), —1 ii p = 2 (mod 4), and zero otherwise. After a 
httle algebra we get 

^ fcgzw y )_ 

This is identically equal to 1 by Euler's pentagonal identity. 



3.2 n = 1, dilute phase. 

This should correspond to the unitary CFT which describes the scaling limit of the critical 
Ising model. Now x = ~^/3, g = 4/3. The numerator in (^ is 



smg 



Now 



p = Qk gives q^'^'' ~^'' 

12A;2-7fc+l 



12fc2+5A:+i 



p = 6k-2 

p = 6k + 1 : q 

p = 6k + 3 

Using the Rocha-Caridi character formula ^H] 

xM = 11(1 - ,0-' E (.^^^^^S^ - {« - -}) , 

r feeZ ^ ^ 

for the case c = |, we see that the first 2 terms give Xi,i and the second pair give Xi,3- 
This then agrees with the result P] for the Ising model with free boundary conditions 

Z = Xi,i + Xi,3 • 

Alternatively we can look at the dual spins, which are fixed on the boundary. If they are 
fixed into the same state on both boundaries we must have p even, so that Z = xi i, and 
if they are fixed into opposite states p must be odd, so Z = xi 3. These also agree with 
Ref. PJ. 



3.3 n = 2 

In this case x = ^iid g = 1, so sin(p + l)x/ sin x — ^ P+ 1- The numerator in (P) becomes 

2 2 

^{p + l)q'^ =^g'r. 



This agrees with the interpretation as the XY model at the Kosterlitz-Thouless transition: 
the terms with p ^ correspond to a total vorticity ±p along the annulus. 
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3.4 Q = 3 Potts model. 

Now X = ^/6 and g = \- The numerator in (^ is 

Esin((»+ l)7r/6) sp^ p 
— ^^^ — -hr^ <!'' '' ■ 
smTT/o 
p 

If we take free boundary conditions on both boundaries we should restrict p to be even. 
Then 





p= 12k 


■Mk^~k 


p 


= 12k + 2 


2„30fe2+9fc+| 


V 


= l2k + A 


30fc2+19fc+3 


p 


= l2k + Q 


_ 30k^+29k+7 


p 


= l2k-A 


_2^30fe2-21fc+/racll3 


p 


= 12k-2 


_ 30A:2-llfe+l 



These pair up as follows: ((1, 6), (2, 5), (3, 4)) to give 

Z = Xi,i + 2X1,3 + Xi,5 , 

which agrees with Ref. jU]. 

Note that if we choose free boundary conditions on one edge and fixed on the other, p is 
restricted to be odd, and the leading term as g — > comes from p = 1, and is 

Z~ V3q^ . 

The v^ is to be expected, because in the Fortuin-Kasteleyn representation each closed 
loop carries a factor a/Q- 



n = Q = 1, dense phase 

In this case x = '^/'^j 9 = 2/3. The numerator in (^J is 

Esin((n+ l)7r/3) Z_£ 
1^ Q 6 6 
sin(7r/3) 

We get a non-zero contribution in the following cases: 



p = 6r 


^br'-r 


p = 6r — 2 


_ 6r2-5r+l 


p = 6r + l 


^6.^+. 


p = 6r + 3 


_ 6r2+5r+l 
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Using Euler's identity again we see that Z = 2, consistent with the dual interpretation 
as the Ising model at zero temperature. (The factor 2 is due to the global spin reversal.) 
On the other hand this model can be interpreted as the Q = 1 Potts model (percolation). 
Choosing the sites on both boundaries to be in the same Potts state enforces p to be even, 
and then we get Z = 1 as expected. 

n = Q = 0, dense phase 

Now X = 7r/2, g = \- The numerator in (^ is 

2_] sin(p + l)vr/2 q~» * 
p 

so p is even. For 

p = 4k : g2fc2-fc 



2 
_ F 
4 



p = Ak + 2 : -q 



2k^+k 

5 



SO 



k 

This is correct, since in this case there is just one macroscopic loop (or spanning tree) 
which is counted with weight n = 0. 



4 Some new results. 

4.1 Percolation. 

By setting cos x = in (^ with (? = | we suppress all other contributions with a non-zero 
number of loops wrapping around the annulus. In terms of percolation, this happens if 
and only if there exists a cluster connecting the two boundaries. This crossing probability 
is therefore 

°° 2 

r=l p 

oo 

r=l feeZ 

so that 1 — P ~ q^^^ as g -^ 0. Using the Jacobi triple product formula this can be written 
in terms of the Dedekind function //(r) = q^^"^^ YY^ii^ ~ q^) with q = e"^'^*/^ as 

^ ^ r^(-l/3r)r^(-4/3r) ^ /, r^(3r)r/(3r/4) 

r/(-l/r)r7(-2/3r/) ^ ^ ' ^(r)r7(3r//2) " 
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In the opposite limit, using ()11|) or the above, we have 

P~ (3/2)1/2 g*, 

as g — > 0, which is the well-known 'magnetic' exponent |T| for the Q = 1 Potts model 
(also known as the 1-arm exponent [201 in the SLE literature.) Note that this result is 
different from, and much larger than, the result found in Ref. JTj. This is because in that 
paper crossing clusters which also wrap around the annulus were disallowed. It would be 
interesting to compare the above result with the implicit formula derived by Dubedat [21] 
using SLE methods. 

Note that, in principle, one can solve for e*^ as a function of n' and substitute in ([T]), hence 
obtaining the complete generating function for the probabilities that a given number of 
clusters wrap around the annulus. 



4.2 Self-avoiding loop: dilute case. 

If we take the 0{n') term in ([Q) with (7 = | we obtain the partition function Zi for a 
single self-avoiding loop which wraps around the annulus. From (Q) we need 



d sin((p + l)x') 



dn' sin x' 



= -i(p+l)cos((p + l)7r/2) 



So we get a non-zero result only when p is odd, say p = 2k — 1, whence, after a little 
algebra. 



2 



Z, = JJ(1 - g')-i ^ A; G ZA;(-l)'=-i g^-^+s . (H) 

r=l k 

The leading behaviour as g — > comes from k = 1 and is 

5 
Zi ~ gs , 

as expected. 

In the opposite limit we can use ()10|). In this case the leading behaviour comes from 
differentiating the exponent: 



2 smx gir'^ ovr 



In g = — I In g I 



If the annulus is mapped into the region between two circles radii ri and r2 > ri, the last 
factor is just ln(r2/ri). 
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4.3 Self-avoiding loop: dense phase. 

In this case 

^1 = -^-'1'' 11(1 - gO"^ E ^ ^°<(P + 1)^/2) fi-"' ■ 

V 

If we let p ^ 2 — p we get the same expression except (p + 1) ^ (3 — p). So the sum is 

- E cos((p + l)7r/2) gV-f , 
p 

and finally (since c = — 2) 

Zi = gi/12 J](l - g^-i E (^'"'"^ - g'''-''+^) . 

r A; 

Using the Jacobi triple product formula this can be rewritten as 

Zi = g-^ J](l-g— 5)2. 

m=l 

The leading term as g — > is 

Z\ ~ g~24 . 

This is reasonable since the loop is weighted by a factor /i^^'^'^s*'^ so the contribution grows 
exponentially with i. 

4.4 Logarithmic cases. 

Logarithmic CFTs have been studied for some time, although the question of how they 
satisfy modular invariance has not been resolved in general ^2- The limit n — > of the 
Q{n) model affords an example in both the dilute and dense regimes; other examples have 
been discussed in Ref. ^2]- If we differentiate the whole expression for Z wrt n at n = 0, 
we get 3 kinds of contribution: the first comes from differentiating wrt %'• this gives the 
partition function for loops wrapping round the annulus as found above, and is a regular 
series in q. The second comes from differentiating q'^l"^^ and gives — (c'(0)/24) Ing times 
the usual partition function at n = 0. The third comes from differentiating the exponents: 
using 

d_,gf_ {g- l)p . ^I^.P 
dg^ 4 2 ' 4 2' 

the result, in the dilute case when g = ^, is proportional to 



(Ing) E(P^ + 2p) sin((p + l)7r/2) q" 
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As before, the only contributions come from p even, and proceeding as before we get 
Inq 5]](A;(2A; + l)q'^^'+^ - k{2k - i)g6fc2-5fc+i) ^ 

k 

or 

k 

Note that we still get the null state structure in this logarithmic sector. 

In the dense phase, the partition function vanishes. The non-logarithmic terms have 
already been evaluated. The logarithmic term is very similar to the above: 



or 



Ing Y.^k{2k + l)q^'''+'' - k(2k - l)q 

k 

Ziog oclng J] k{2k + l)(g2'='+'= _ q 



2fc2-3fc+l^ 



2k^+k 2fc2+3A:+l^ 

k 

once again showing the null states. The leading term as g -^ gives the contribution of 
single dense loops (or spanning trees) which do not wrap around the annulus: note that 
this is much smaller than the 0{q~^f'^'^) contribution from those which do. 

5 Summary and further remarks. 

In this paper we have presented an explicit result for scaling limit of the partition function 
of the critical 0(n) and Q-state Potts models on the annulus. Our formalism makes it 
simple to count loops which wrap around the annulus with different weights, leading 
potentially to many new formula for crossing probabilities in percolation and for self- 
avoiding loops, some of which have been presented here. 

The electric-magnetic dual of our arguments for a long (£ ^ L) annulus can in fact 
be applied to a long (L ^ I) cylinder, to give an alternative derivation of the usual 
relation between g, x ^"^^ ^- For a long annulus we argued that ±2 magnetic charges 
(vortices) would rearrange themselves in such a way as to induce mean magnetic charges 
d=mo = ±(1 — g)/ g at the ends of a long rectangle, leading to a net magnetic flux around 
the annulus. For a long cylinder, one may similarly argue that the ±2 electric charges in 
the model also rearrange themselves to give net electric charges ±eo = ±(1 — g) at either 
end. These can then be interpreted as counting loops going around the cylinder with the 
weight n = 2cos7reo. 

Like all Coulomb gas methods, however, the arguments are somewhat heuristic, although 
they lead to completely explicit formulae, and, because of the complex weights, it seems 
hard to make them rigorous. In particular, although we have computed the partition 



function, it is by no means clear that the same ensemble can be used to compute correla- 
tion functions in the original model. It would be nice to see a direct connection with the 
other 'Coulomb gas' approach which has been employed in CFT, namely that originally 
developed by Dotsenko and Fateev [SH!- This is essentially a way of constructing holo- 
morphic conformal blocks using modified free field theory. However, it has many features 
in common with the Coulomb gas construction used here: the background charge mo, 
and the marginal screening operators, which in our case are the ±2 vortices. In the bulk, 
it is still necessary, in the Dotsenko-Fateev approach, to sew together the holomorphic 
and anti-holomorphic blocks in a consistent way to obtain correlation functions, but in 
boundary CFT the correlation functions are linear combinations of the conformal blocks 
(specialised to real values of their arguments), and so the correpondence between the two 
approaches should be more direct. It would, of course, be important to establish any of 
these results rigorously, for example by using SLE methods. 

Acknowledgement. This work was supported in part by EPSRC Grant GR/R83712/01. I 
thank H. Saleur for pointing out several important references missing in an earlier version 
of this paper. 



A Boundary terms in the gaussian model. 

Consider the action 

S = So + Si = — / {yKfdxdy + ai dyh{x, y = 0)dx + 0:2 / dyh{x, y = L)dx . 

We wish to compute the regularised free energy, or equivalently the ground state energy 
Eq of the associated hamiltonian. Let 






S,= j dxY^{\fl + \{n7^/Lffl^ 

n 



Then 

from which we read off the ground state energy 

Eo = -^u{nn/L), 

n 

where to'(fc) = k. We can regularise this sum either by modifying the dispersion relation (eg 
using a lattice, in which case u;(/c) = 2sin(fc/2)) and using the Euler-Maclaurin formula, 
or using zeta-function, in which case we get the standard result 

So = ^C(-l) ' 



2L'' ' 24L 
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corresponding to c = 1. 

Now add in 

5i = (47rM)^/2 J](n7r/L)(ai + (-l)"a2)/n. 

n 

Completing the square, the contribution of the nth mode can be written 

so the change in the ground state energy is 

El = -{2tx /gL) ( 5^(«i - tta)' + ^ («i + as)' j • 

\nodd neven / 

If we use zeta function regularisation we have 

Y^l = hm(5^n--5^(2n)-) = hm(l-2-)C(s)=0 

nodd n n 

El = li-B2-r = c(o) = -^, 

neven n 

SO 

Ei = ^{ai + a2f. (12) 

If we use a lattice dispersion relation and also replace the boundary derivatives by finite 
differences, the sum becomes 

^ cos2(n7r/2L) ' 

n=l ^ ' ' 

on which we use the formulae 



5^/(n/L) = \l\ f{x)dx + 0{L- 



1-1 m rfiTi '-' 



/(^)c?a;-— /(0) + O(L-^) 

neven 



giving the same result. It can also be verified by writing, for a general position dependent 

a 

Z = Zo /exp ( f a{l)d^h{l)dl\ \ = Zq exp (- [ [ a{l)a{l')d^d'^G{l, l')dldA , 

where G is the Green's function for the free field with Dirichlet boundary conditions. 
fjl2|) leads to a modification to the effective central charge 



24 

c = 1 (ai + 02)', 

9 



as claimed in Sec. El 
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